Many-body effects in magnetic inelastic electron tunneling spectroscopy 
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Magnetic inelastic electron tunneling spectroscopy (lETS) shows sharp increases in conductance 
when a new conductance channel associated to a change in magnetic structure is open. Typically, the 
magnetic moment carried by an adsorbate can be changed by collision with a tunneling electron; in 
this process the spin of the electron can flip or not. A previous one-electron theory [Phys. Rev. Lett. 
103, 176601 (2009)] successfully explained both the conductance thresholds and the magnitude of the 
conductance variation. The elastic spin flip of conduction electrons by a magnetic impurity leads to 
the well known Kondo effect. In the present work, we compare the theoretical predictions for inelastic 
magnetic tunneling obtained with a one-electron approach and with a many-body theory including 
Kondo-like phenomena. We apply our theories to a singlet-triplet transition model system that 
contains most of the characteristics revealed in magnetic lETS. We use two self-consistent treatments 
(non-crossing approximation and self-consistent ladder approximation). We show that, although the 
one-electron limit is properly recovered, new intrinsic many-body features appear. In particular, 
sharp peaks appear close to the inelastic thresholds; these are not localized exactly at thresholds 
and could influence the determination of magnetic structures from lETS experiments. Analysis of the 
evolution with temperature reveals that these many-body features involve an energy scale different 
from that of the usual Kondo peaks. Indeed, the many-body features perdure at temperatures much 
larger than the one given by the Kondo energy scale of the system. 

PACS numbers: 68.37.Ef,72.15.Qm,72.10.Fk,73.20.Hb 
Keywords: lETS, spin-flip, Kondo effect, singlet-triplet Kondo 



I. INTRODUCTION 

Magnetic inelastic electron tunneling spectroscopy 
(magnetic lETS) detects magnetic excitations on a sur- 
face by measuring the changes of conductance of a scan- 
ning tunneling micrsocope (STM) junction when an ex- 
citation is producedi. As in general lETS, the possibility 
to excite the surface leads to an increase of the number 
of possible final channels for the tunneling electron when 
the junction bias matches an excitation energy thresh- 
old, consequently leading to an abrupt increase of the 
junction conductance. This high sensitivity on the small 
magnetic energy scale has permitted Hirjibehedin and 
coworkers^ to measure the magnetic anisotropy energy 
(MAE) of single atomic adsorbates as well as the mag- 
netic coupling among thent^. More recent experiments 
revealed the change of sign of MAE when a phthalocya- 
nine molecule is adsorbed on a surface^ and the existence 
of an exchange coupling between magnetic molecules in 
a multilayer setup*^. This technique has open a venue 
to the study of magnetic phenomena at surfaces on the 
single atom or molecule scale. 

Theories have been developed to rationalize the steps 



in conductance found in magnetic JETS with a good de- 
gree of success^"—. The large changes in conductance ob- 
served in magnetic JETS were explained by the large par- 
ent coefficients of the initial and final adsorbate states in 
the tunneling state — . During tunneling the spins of the 
tunneling electron and of the magnetic adsorbate cou- 
ple together into a total spin St that characterizes the 
spin symmetry of the tunneling process. The excitation 
process is then pictured as the coupling/decoupling of 
the two spins. Since the parentage coefficients of the St 
states in the initial and final magnetic states of the adsor- 
bate can be large (these are Clebsch-Gordan coefficients 
associated to MAE structure coefficients), the probabil- 
ity of forming various excited final states can be very 
large. In this way, the flux of incident electrons is shared 
among the possible final channels that are energetically 
possible; the branching ratios among channels are sim- 
ply governed by structure coefficients. This theory also 
showed that the magnetic excitation of an adsorbate may 
or may not imply spin-flip of the tunneling electron^. 

The Kondo effect is associated with a phenomenon 
bearing many links with magnetic JETS: the fluctuations 
induced by the spin-flip of an electron during its collision 
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with a magnetic impurity at constant energy. Kondo 
physics and inelastic effects have been extensively stud- 
ied in quantum dots and nanotubesi^ri^. Recently, vi- 
brational side-bands have been predicted and reported 
for molecules displaying Kondo peaks at zero biaai^"— . 
Inelastic processes in Kondo physics have been studied by 
Zarand, Borda and coworker a^^i^^ who have shown that 
the usual Kondo theories can reveal the amount of elastic 
and inelastic spin-flip in the electron-impurity scattering 
event. For the particular case of magnetic lETS, Zitko 
and Pruschke^^ have applied Kondo theories to the study 
of the coexistence of a Kondo peak and of lETS steps 
in the STM conductance of Co atoms on CuN/Cu(100) 
substrates. Unfortunately, the numerical renormaliza- 
tion group method that was used is not very accurate 
in reproducing the sharp conductance steps. However, 
these authors gave the first unified picture of Kondo and 
magnetic lETS. 

More recently. Hurley and collaborators^! have used 
the Kondo Hamiltonian and perturbation theory to ana- 
lyze the magnetic lETS of Co and Fe on CuN/Cu(100). 
The authors conclude that certain spike- like structures at 
the inelastic thresholds in experimental lETS are actu- 
ally due to a Kondo-like effect. However, spike-like struc- 
tures close to the inelastic thresholds can also be found 
due to electronic heating under high current conditions 
as recent experimental and theoretical rcports^^"— show. 
Indeed, the conductance changes in a non-linear way if a 
tunneling electron probes the adsorbate still excited by 
earlier tunneling electrons. 

Here, we study the magnetic transitions between sin- 
glet and triplet configurations of a magnetic adsorbate. 
We use our previous study on adsorbed copper phthalo- 
cyanine^^ and parametrize it to study the role of different 
ingredients in the spectral function. Our study considers 
two independent spins localized in two different molecu- 
lar orbitals whose interaction is simply given by a Heisen- 
berg term of exchange interaction /. This model contains 
similar physics to the recent two-impurity Kondo system 
studied by Bork and co-workers2^. 



II. ONE-ELECTRON THEORY 

The magnetic excitation of magnetic impurities on a 
solid surface using an STM was modelled as an electron- 



impurity collision in Refs. [9l l27ll31 |. During the colli 



sion of the electron with the impurity, the tunneling 
electron and adsorbate spins briefly couple together. A 
transient resonant state with a finite lifetime can be in- 
volved, but not necessarily; it could simply be a scat- 
tering state that fixes the tunneling symmetries. These 
symmetries depend on both the electron energy and the 
electronic structure of the impurity. In the present case, 
the electron energy is very small since typical magnetic 
excitations are in the meV range. Hence, the impurity's 
electronic structure at the Fermi energy is the relevant 
one during the collision. In Refs. [9l l27H3]| the impurities 
had orbitals straddling the Fermi energy of the substrate 
which were thus involved in the tunneling process. 

In the present work, we extend the study to systems 
where a positive ion is a more likely origin of the transient 
state. In particular, we are studying the singlet-triplet 
excitations of a magnetic impurity with a large charging 
energy U. In this case, the negative ion is energetically 
less accessible than the positive one. This feature has 
been used to model the Kondo effect in certain molecular 
system o^^i'^^" — . 

The modelling of the Kondo effect runs parallel to the 
one-electron theory used to account for magnetic excita- 
tions in impurities. Indeed, the Kondo effect (see for ex- 
ample Ref. [28j and references therein) can be described 
as a fluctuation between a charged transient state and 
the ground state of the impurity. The Kondo effect builds 
up on the coherence between impinging electrons, it is a 
genuine many-body effect, while our earlier lETS stud- 
^Qg9 ,27,3i pjjjy consider single-electron collisions. 

Both in the one-electron and many-body pictures, the 
tunneling state connects the initial state of the impurity 
with a final state that can be different, leading to mag- 
netic excitations. We study below how many-body effects 
influence magnetic excitation processes. 



Inelastic transition rate 



In the first part of this work, we use our previous 
one-electron theorj*^ to evaluate the magnetic lETS for 
a singlet-triplet excitation. In the second part, we use 
the non-crossing approximation^^ and the self-consistent 
ladder approximation^ for the same model system. The 
sharp behavior of the conductance steps due to lETS is 
accounted for by both theories. However, many-body ef- 
fects appear to be very important at the inelastic thresh- 
olds. Depending on the parameters of the impurity, we 
further find that the actual excitation thresholds can 
be substantially shifted when many-body effects are in- 
cluded. This finding can have important consequences in 
the determination of MAE and more generally of adsor- 
bate magnetic structures from lETS. 



We assume the tunneling process (electron-adsorbate 
collision) to be very fast, much faster than the interac- 
tion at play in the singlet-triplet splitting, so that we can 
resort to the sudden approximation. The T-matri»2Zi^ 
between an initial state \i) and a final state |/) of the 
complete electron-l- adsorbate system is obtained from the 
sudden T-matrix expressed in the basis set of the inter- 
mediate states \St, Mt) (spherical spin symmetry is kept 
during the brief collision). In the present case, since the 
system ground state is a singlet, the tunneling symmetry 
can only be S't=1/2. 

In the large U limit, the collision takes place between a 
hole and the impurity. The ground state of the molecule 
is a singlet, then the initial state is a singlet times (ten- 



serial product) a hole of the conduction band. Hence, 

K) = c.|o,o) (1) 

where Ca destroys an electron with spin cr, and the impu- 
rity in the ground state is given by its spin |5, M) which 
is a singlet |0, 0) . Similarly, the final state is 
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\S,M) 



(2) 



where the impurity is left in one of the jS", M) states and 
the hole is in a a' state. Below, we assume the STM tip 
to be unpolarized, so that we sum over the contributions 
for cr and a' spins. 

The sudden T-matrix reduces to a projection operator 
on all the \St, Mt) states times a common transmission 
amplitude, B^'^. In tunneling, the transmission proba- 
bility density is proportional to the density of states of 
the sample, p{uj). This still holds in more complicated 
situations as shown in Ref. [s^. Then, the transmission 
probability density is proportional to a constant coming 
from the transmission amplitude, |-Bp, times the density 
of states. Hence, the elastic probability density for an 
electron energy ui, Tg, is given by 

T,H = y: I E ^H5.io, o)n^ = 

cr Mt 

(3) 

The inelastic contribution, Ti, contains three possible 
orientations, Mf, of the impurity's spin because it is a 
triplet. Hence, 
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a,cr',Mf Mt 



Y,B{QMcI\St.Mt) 



X {ST,MT\c„,\l,Mf)\^ 



(4) 



At this point, we can stress that the inelastic conduc- 
tance is associated to both spin- flip (ct ^ a') and non 
spin- flip ((T = cr') processes for the tunneling electron; ex- 
citation of the triplet state by a tunneling electron with- 
out a change of the electron spin is 50 % less probable 
than excitation of a triplet adsorbate state with a spin- 
flip of the tunneling electron. 

From Eqs. ([3]) and Q the relative contribution of the 
elastic part to the transmision probability density at an 
inelastic threshold is just 



(5) 



and 3/4 is the relative contribution of the inelastic part 
of the transmission probability density. 

In the spirit of the generalization of the Landauer 
transmission formula to include inelastic transitions^^—, 
it is possible to link2^ the transmission probability den- 
sity to the projected density of states on the magnetic 
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FIG. 1: Projected density of states on the magnetic adsorbate 
electronic structure or spectral function A as a function of the 
electron energy u) with respect to the Fermi energy of the sub- 
strate. The vertical dashed lines show the one-electron inelas- 
tic thresholds. The figure presents the adsorbate density of 
states in the absence of inelastic efi'ects (marked as "No inelas- 
tic"), the many-body spectral function and the one-electron 
conductance with inelastic effects. The system's temperature 
is 7 K. 



adsorbate or spectral function. A, and we use, below, this 
link to compare the results of the one-electron and many- 
body studies. In the present one-electron approach, the 
spectral function. A, is obtained as: 

A = p{^){-^ + -^[f{u + I) + f{I^uj)]), (6) 

where / is the singlet-triplet excitation energy and f{uj) 
is the Fermi function of the substrate. 

Fig. [T] shows the characteristic step-like function for 
the spectral function, A, at low temperature. The two 
steps at positive and negative energy are associated to 
the opening of the inelastic triplet channel. The results 
in Fig. [T]use the parameterization of the adsorbed CuPc 
molecule determined in Ref. [28l . in particular the excita- 
tion energy, /, is equal to 25 meV. 



B. Thermal effects 

The above one-electron results in Fig. [1] correspond to 
a low temperature of the surface. Thermal effects tend to 
wash out the inelastic structures via two contributions: i) 
the Fermi function is not a step function and this rounds 
the conductance steps at threshold and ii) at equilibrium 
at finite temperature, the system is not initially entirely 
in the singlet state, a small thermal population of the 
triplet is also present. 

For the finite-temperature differential conductance, 
the first effect has been discussed in detail in Refsi '^^'^^ . 
The abrupt steps at inelastic thresholds that are present 
in the conductance at vanishing temperature are replaced 
by rounded steps due to the smearing of the energy distri- 
bution of the electrons at finite temperature. This leads 
to a significant broadening of the step function of the 
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order of 5.5 ksT—. In the present one-electron results 
however, for a consistent comparison with the many body 
results on the spectral function, we evaluate the above 
one-electron spectral function, i.e. we do not include the 
broadening effect coming from the thermal distribution 
of the electrons in the tip. 

The second effect is due to the finite thermal popu- 
lation of the triplet state by the Boltzman factor, F = 
exp[— J/fefiT] and the corresponding decrease of the sin- 
glet population. The total spectral function. A, is then 
equal to: 



the substrate, 



E 

k,acr 



(8b) 



The substrate electronic degrees of freedom are spin 
tr = ± 1 , channel a = 1,2 and the remaining degrees 
are encapsulated in the k symbol. 
The impurity Hamiltonian 



h — e \aa){aa\ + /Si • S2 



(8c) 



A = 



{| + |[/(- + /) + /(/--)]}/"^ 
1 + 3F 



3F 



(7) 



Note the factor 3 due to the triplet degeneracy. The 
total (elastic -|- inelastic) contributions of the singlet and 
triplet states are equal. Since the triplet contribution 
does not exhibit any step, the thermal population of the 
triplet state tends to smooth out the stepped structure. 
This effect is only visible on the results in Fig. [5] at the 
highest temperature. 



III. MANY-BODY THEORY 

The extension of the above one-electron theory to the 
many-body case is achieved by keeping the coherence be- 
tween impinging holes. This is a tremendous task, but 
can be easily achieved using self-consistent schemes such 
as the non-crossing approximation (NCA) and the self- 
consistent ladder approximation (SCLA). 



Non-crossing approximation: Anderson 
Hamiltonian 



We consider that the impurity is fluctuating between 
two charged states: one corresponding to the one-body 
ground state, and one of the above transient states. The 
source of fluctuation is the hybridization of the impurity 
with the substrate, V. By allowing the two charge states 
to evolve self-consistently, NCA is an all-orders theory, 
albeit neglecting all terms that "cross" . It is at the sixth 
order in V that the first crossing terms of the pertur- 
bation expansion are neglected. NCA is then a method 
for the solution of the Anderson Hamiltonian, where spin 
fluctuations are brought about by the hybridization term, 
V. 

The considered Anderson-like Hamiltonian contains 
three terms, 

H = Hs + V + h. (8a) 
The first term is the free-electron-likc Hamiltonian of 



has been used to describe a metal-organic adsorbate^^. 
The doubly degenerate ligand orbitals have orbital in- 
dex a = 1, 2 and are represented by Hubbard operators, 
which project out configurations where ligand occupation 
is higher than one. The ligands are subject to exchange 
interaction with a third orbital strongly localized at the 
molecular center which is represented by a spin-half oper- 
ator 82- Its charge does not fluctuate because it is a very 
compact orbital decoupled from the metallic substrate^^. 
The ligand spin operator can be expressed through the 
vector of Pauli matrices r as follows 



a=1.2 



The substrate - impurity hybridization is expressed by 
the term 



^ = E (^k4,aJ0)(aa| + Vk*|aa> 



(8d) 



k.ao" 



which docs not mix spin, a, and orbital ,a, degrees of 
freedom. The empty ligand configuration is denoted by 
|0). 

Hence, Hamiltonian ([8a|) describes^^: (z) an electron in 
an orbital disconnected from the reservoir, (ii) the charge 
fluctuations of the 2-fold degenerate orbital connected to 
the reservoir, and (Hi) the mutual interaction between 
both electrons via an exchange term of matrix element 
/ > 0. This system can have singlet-triplet excitations 
with an excitation energy equal to I. 

The substrate-impurity hybridization enters in the 
NCA equations via the energy-dependent width of the 
impurity orbitals. In the present work, we model it by 
a rectangular function that is zero for electron energies 
beyond the electron band, of bandwidth 2D. 

For / = the spin S2 decouples and we are left with a 
SU(4) Kondo effect with a Kondo temperature given 
bj^: 



(9) 



where 27rr is the level width of the orbital a resonant 
with the metal substrate. sets a natural energy scale 
of the present problem and we will use it as the energy 
unit of the different calculated quantities. 
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B. Self-consistent ladder approximation: 
Coqblin-Schrieffer Hamiltonian 

It is also interesting to consider the above physical 
model adapted to the Kondo Hamiltonian. The Kondo 
Hamiltonian explicitly includes a spin-spin interaction 
term between an itinerant electron and the impurity spin. 
In the Kondo limit, |e| ^ F, where the impurity or- 
bital energy, |e|, is much larger than the level broad- 
ening, r due to the hybridization term, V, both An- 
derson and Kondo Hamiltonian describe the same spin- 
flip physics^i^. The Coqblin-Schrieffer Hamiltonia n^^ 
generalizes the Kondo Hamiltonian to include orbital de- 
grees of freedom, generalizing the SU(2) Kondo problem 
to SU(N), where N is the combined orbital and spin de- 
grees of freedom of the impurity. The Coqblin-Schrieffer 
Hamiltonian is interesting to be considered here because 
it allows us to explore spin excitations in the Kondo 
limit. We will use the equivalent of NCA for the Coqblin- 
Schrieffer Hamiltonian, namely the SCLA^i^. 

The Coqblin-Schrieffer Hamiltonian can be obtained 
from the Anderson Hamiltonian, Eq. (|8]), in the 
Schricffer- Wolff limit 



-oo, — = const. 



(10) 



Hence, Hamiltonian (|5]) can be approximated by the 
Coqblin-Schrieffer Hamiltonian^: 



(11) 
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C. Equilibrium regime 

The following results have been obtained using equi- 
librium NCA and SCLA. In principle, a non-equilibrium 
calculation is needed to account for the correct coherence- 
decoherence balance at the excitation bias^. However, 
the situation studied here is the one found in STM studies 
of molecules on surfaces. Typical parameters are tunnel- 
ing currents below the nA range and bias of a few meV. 
Assuming a single impurity resonance, and a small bias, 
Vtip, such that the current takes place in resonance, the 
current can be estimated by a Breit-Wigner-like expres- 



r-^r 



tip 



(12) 



where Ytip is the resonance broadening due to hybridiza- 
tion with the STM tip. If we take typical STM param- 
eters such as 1 nA, 0.1 V, we obtain from Eq. p2)) that 
r^jp ss 10/xey. Typical F for molecules on surfaces are 
larger than 100 meV (see for example the calculations 
of Ref. ) . The non-equilibrium modification of NCA 



equations^ for slowly varying substrate density of states, 
reduces to the replacement of Fermi occupation func- 
tions, f{uj — fi), by an effective distribution function Feff 
given by^: 



tip 



Ttip + r 



/(t^-Mtip)' 



-J{LU-fl). (13) 



Since F^^p is easily a factor 1000 smaller than F (a factor 
10^ in the above example), we recover -Fe// ~ /("^ ~ A')- 
Hence, in typical STM inelastic measurements, the sam- 
ple will be largely in equilibrium. Qualitatively, the tip is 
extracting only an extremely small electron current from 
the molecule, much smaller than the electron fluxes that 
come from the substrate electron bath. If one further 
note that the bias applied to the junction is very small 
(tens of meV) compared to the local electrostatic poten- 
tials, then one can assume that the electronic structure 
of the probed molecule is not modified by the presence 
of the tip. This is the customary Tersoff-Hamann pic- 
ture^i^, where the STM is assumed to read the unper- 
turbed spectral function of the molecule on the substrate. 



D. Singlet-triplet excitations 

Reference [1^ gives a detail account for the implemen- 
tation of NCA in the case of singlet and triplet molecules. 
As we showed in the previous section, the Hamiltonians 
contain an exchange term between two spins localized in 
the molecule: 



Hj — /Si ■ S2 



(14) 



The two spins relate to two different molecular orbitals 
that couple differently with the metallic continuum. This 
gives rise to a rich variety of physical situations depend- 
ing on the value of the exchange inter ation /, and the 
Kondo scales of the different orbitals (see discussion in 
Ref. [2^ and references therein). 

Here, we just consider positive /, such that the molecu- 
lar ground state is a singlet and the singlet-triplet energy 
excitation is equal to /. The case of / < exhibits fea- 
tures similar to the ones described in the present paper 
(see Ref. [2^ ) . Figure [1] shows the results of the singlet 
molecule for an exchange interaction / = 25 meV— ; the 
Kondo temperature, T]^, is equal to 30 K. These results 
show that the general features of the inelastic effect, i.e. 
the sharp conductance steps at inelastic thresholds, are 
readily understood with the above one-electron theory. 

However, new features appear at the excitation thresh- 
olds that are purely of many-body character. Mainly, the 
spectral weight is greatly increased near the thresholds. 
This is reminiscent of the Kondo peak at zero energy. The 
difference is that in a SU(N) Kondo peak, there are fluc- 
tuations between degenerate impurity orbitals induced 
by spin-flip transitions at constant energy; in contrast, 
here, the singlet-triplet fluctuations involve an energy 
change of the adsorbatc that is provided by the junction 



bias. In addition, in the present case, the singlet-triplet 
transitions are not pure spin-flip transitions, they can 
also occur without a change of the spin of the electron 
colliding on the impurity. 

The structure appearing at zero energy in the many- 
body spectral function is due to the self-consistent ap- 
proaches used here. The reason for its appearance is the 
artificial fiow of the marginal potential scattering term in 
NCA4^i^, which overestimates potential scattering and 
hence leads to a spurious Kondo-like feature at zero en- 
ergy. 

The results in Fig. [T] were obtained using the 
parametrization from Ref. [l^ for the adsorbed CuPc. 
From now on in this paper, we will vary the parameters 
in this modelling (ratio of excitation energy and Kondo 
temperature, energy of the orbital) in order to decipher 
the role of the various parameters in the characteristics of 
the many-body features appearing close to the inelastic 
thresholds. 



E. Temperature effects 

When the temperature, T, is much smaller than the 
excitation energy (fc^T ^ /) thermal effects due to the 
equilibrium population of the triplet state are negligible. 
The one-clcctron cases of Fig. [T] and to a lesser extent of 
Fig. [2] are in this regime. Hence, the dominant thermal 
effect is due to the smearing of the sharp conductance 
steps at the inelastic thresholds due to the Fermi function 
broadening. 

Figure [2] is a systematic study of the temperature ef- 
fect in the spectral function with and without many-body 
effects for / = 10 fcsT]^. As for Fig.[Tl the many-body re- 
sults have been obtained for the Anderson Hamiltonian, 
Eq. ((5a)) solved using NCA. The many-body peaks appear 
to vary rapidly with T, they collapse as the temperature 
is increased. This behavior, for moderate temperatures 
(kBT <^ I), is due to the excitation of thermal electron- 
hole pairs that destroy the electron coherence, reducing 
the rate of coherent singlet-triplet transitions and hence 
the enhanced density of states at the excitation thresh- 
olds, similarly to the disappearance of zero-energy Kondo 
peaks with temperature. 

The comparison of the one-electron and many-body 
spectral functions shows that the apparent energy thresh- 
olds are displaced one with respect to the other. This 
is already apparent in Fig. [1] where the threshold coin- 
cides with the mid-point of the inelastic step for the one- 
electron spectral function but the mid-points of the in- 
elastic steps in the many-body results are clearly shifted 
from the energy threshold. Figure [2] shows a common 
point where the one-electron spectral functions cross for 
the five considered temperatures. However, there is no 
common crossing point in the many-body curves and 
their behavior is controlled by the temperature evolution 
of the Kondo-like peaks. 

For the SU(4) Kondo effect (case of a vanishing ex- 



0.006 




FIG. 2: Temperature dependence of the spectral function for 
an exchange coupling of the impurity's spins equal to ten 
times the Kondo temperature, / = 10 ksT^, for the An- 
derson model, Eq. (|8a|) . The energy axis is expressed in 
units, such that the inelastic thresholds appear at ±10. As 
the temperature is raised (from 0.2 to 2.4 T^), the Kondo- 
like peaks located close to the inelastic thresholds are smeared 
and finally even the inelastic gap disappears. Dashed lines are 
the results for the one-electron theory, that accounts for the 
inelastic changes of the spectral function where all impurities 
are assumed to initially be in their singlet state. 



change interaction, /), the Kondo peak is greatly dimin- 
ished at T « T^. This is seen in Fig. [3] (a), where 
the / = 0, SU(4) Kondo peak is plotted for several 
temperatures. Figure [3] has been obtained by solving 
the T-matrix, T, for the Coqblin-Schrieffer Hamiltonian, 
Eq. (|TT|) . The figure shows the imaginary part of the 
T-matrix times the density of states p as a function of 
the electron energy w. This is equivalent to plotting the 
spectral function as a function of w, since the hybridiza- 
tion function F is the proportionality factor connecting 
them. Comparison of Figs. [2] and [3] shows the equiva- 
lence of both Hamiltonians, Eqs. (|5a)) and ([TT|) . and of 
the solution methods. 

Surprisingly, when the temperature behavior of the in- 
elastic Kondo-like peaks is studied (Fig. [3] (6)), we ob- 
serve that the peaks survive the temperature increase 
much longer than the SU(4) peaks (Fig. |3] (a)). Figure[3] 
(b) shows —ImpT as in (a), for / — 256 ksT^^. However, 
we see that the spectral features are still important at 
T w T^, and even at T = 50 spectral peaks are still 
visible. However, at T = 2.0 the SU(4) peaks are 
very diminished. Fig. [3Ja). This shows that the inelastic 
process with an excitation energy equal to / sets in a new 
energy scale that controls the spectral variation with tem- 
perature. Only when is very small (results not shown 
here), we recover the one-electron limit. Hence, the con- 
ditions where many-body effects are not observable, while 
inelastic effects are observable are kb -^kbT <^ I. 

For temperatures above J, the two inelastic peaks co- 
alesce giving rise to a a single broad peak, similarly to 
the SU(4) limit. This variation appears in more detail in 
Fig. m which presents the heights of the inelastic Kondo 
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FIG. 3: Imaginary part of the T-matrix times tlie density of 
states p for different temperatures in units, (a) SU(4) 
case for the spin 1/2 problem obtained when the exchange 
interaction between the two localized spins is set to zero, I = 
0. (fe) the singlet-triplet excitation for a large value of the 
exchange integral, I = 256 /cbT)^. 



peaks relative to their value at low T as functions of the 
temperature. 

For 7 = 0, the peak heights are approximately de- 
scribed by the blue curve (Fig. S} given by: 



1 



l + (2iA-l)(T/T0) 



(15) 



with s = 0.4. The same expresion, for s = 0.21 has been 
used to fit the results for the conductance in the Kondo 
regime from a renormalization group calculatio n^^'^*^ . 
Since, the fitted quantities (T-matrix vs conductance) 
and the system (SU(4) vs SU (2)) are not the same, it is 
not surprising that s is different. This fit ensures that the 
peak has dropped to 1/2 at T = T^, which shows that 
the expression for the Kondo temperature, given by 
Eq. ([2]) is a satisfactory approximation for both the NCA 
and SCLA. 

Figure 2] also plots the temperature evolution of the 
heights of the two inelastic Kondo peaks for / = 10 fcsTj^- 
(down triangles and diamonds for lower and uppper 
peaks, resp.) and / = 256 kBT^ (squares and up trian- 
gles for lower and uppper peaks, resp.). Each peak height 
is normalized to its value for T « 0. At T « 6r^, the 
two peaks of the / = lOfc^T]^ case coalesce and form a 
single peak. In order to show this clearly in Fig.|3]we have 
to change the normalization of the peak heights above 6 




FIG. 4: Maximum value of the spectral functions of Figs. [2] 
and [3] as a function of temperature normalized to the value 
at T ~ 0. Black dots: peak near the Fermi energy for I — 
(Fig. El^a)); mauve down triangles: peak near u = — lOfcsT]^ 
and cyan diamonds: peak near cj = -10fcflT£, for I = 
10 ksT^K (Fig-Elli red squares: peak near w = —256 ksT^ and 
green triangles: peak near lu = 256 fcsTj^ for / = 256 ksT^ 
(Fig. [3{fe)). The blue line is an empirical fit to the results for 
1 = that displays at T = a maximum value half of that 
at T = 0. The normalization of the down-triangle curve has 
been changed above T = 6T^, due to the coalescence of the 
two inelastic Kondo- like peaks at high T (see text for details). 



T^. We chose to normalize the single peak above 6 T]^ as 
the high energy peak. This increases the discontinuity in 
the down-triangle curve but preserves the diamond curve. 
Beyond T « 10 the curve evolves more softly. This 
behavior shows that / is the energy scale that governs 
the spectral-feature evolution at high temperatures. For 
the / = 256 ksT^ cases, the plotted temperature range 
(see e.g. in Fig. [3] (6)) does not reach a point at which 
the two peaks overlap significantly. 

The different sensitivity to temperature increase in the 
inelastic case as compared to the 'usual' Kondo problem, 
can be attributed to the different role of electron-hole 
pairs. Decoherence increases as a temperature rise in- 
duces more electron-hole pairs. For the usual Kondo case, 
a Fermi electron is scattered by the impurity, hence ther- 
mal electron-hole pairs are very efficient in causing the 
electron decoherence. However, in the inelastic case, an 
excited electron is scattered by the impurity and so it 
concerns the energy range around I. Hence, decoherence 
becomes particularly efficient when the thermal electron- 
hole pairs have enough energy to reach the inelastic tran- 
sition range, i.e. when the temperature is of the order of 
the excitation energy. This explains why as we increase 
/, the Kondo-like peaks survive at higher temperatures. 



Orbital energy dependence of the excitation 
threhold 



The inelastic spectra show Kondo-like features at the 
excitation thresholds. Zero-energy Kondo peaks are basi- 
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cally determined by the value of the Kondo temperature, 
T^, and of the orbital energy, e. is responsible for 
the zero-temperature width of the Kondo peak and e is 
an important factor for the actual shape of the Kondo 
peak by virtue of the Friedel-Langreth sum rule^i^. For 
this reason, we have studied the variation of the inelas- 
tic features with the orbital energy, keeping the Kondo 
temperature, T]^, constant. 

Figure [S] shows the evolution of the inelastic features 
for five different values of the orbital energy (excitation 
energy / equal to 4 fc^T^ in Fig. ^a) and to 32 ksT^ in 
Fig. El^b)). The first four orbital-energy values increase 
by a factor of two: |e| = 110, 220, 440 and 880 ksT^- 
The fifth value is computed using the Coqblin-SchriefFer 
Hamiltonian, Eq. (fTT|) . that corresponds to the Kondo- 
limit or e — > — oo. The negative-energy peaks of the spec- 
tral functions have been normalized to one by multiplying 
the full spectral function by a constant number. The be- 
havior of the low and high- energy tails can be understood 
just by the change in the density of states as the orbital 
energy shifts, because is kept constant and hence, the 
ratio e/r is constant, where 27rr is the width of the one- 
electron peak originating at e. Hence, as e approaches 
— oo, the width 27rr increases. For |e| — 110 ksT^, 
the orbital resonance is close to the Fermi energy, and 
hence, the density of states drops rapidly, while for a 
larger value, the resonance is far from the Fermi energy 
and much broader, dropping more slowly following the 
trends of Fig. [5j 

Despite the relation of the peaks at threshold with 
Kondo features, it is difficult to conclude on some type of 
extension of the Fermi-Langreth sum rule since the two 
peaks at the inelastic thresholds present different behav- 
iors. 

As |e| increases the peaks at thresholds move away from 
the threshold energy towards positive energies. However, 
the inset shows that the e — > — oo case for / = 32 fc^ is 
anomalous in the sense that it presents an extra broaden- 
ing instead of a positive-energy shift. This hints at some 
saturation effect as the energy scales increase. 

It is noteworthy that the evolution of the negative- 
energy peaks is faster than the positive ones. The peaks 
at negative energy are also narrower than the positive 
energy ones and their relative heights change depending 
on the values of / and e. While for / = 4 fegT]^ the 
positive-energy peak is equal or larger than the negative- 
energy case, for / = 32 fesT]^ they are smaller. 

Figure [5] also shows that the deviation of the peaks 
from the inelastic thresholds. A/, depends on /. We 
study this behavior in the next section. 



G. Asymptotic behavior of the threshold 
renormalization 

All the above results show that the apparent inelastic 
thresholds shift as / increases. Actually, the threshold 
shifts lead to a reduction of the inelastic gap. We can 




FIG. 5; Normalized spectral function as a function of the 
electron energy in units of the Kondo temperature, T^. Case 
(a) corresponds to an excitation energy I — AksT^ as de- 
picted by the vertical dashed line. Case (b) corresponds to 
the / = 32 kBT^ case. The curves correspond to orbital ener- 
gies |e| = 110, 220,440, 800 /cbT)^ for the curves in increasing 
spectral function in the positive-energy part of the graphs. 
The topmost curve has been computed using the equivalent 
Coqblin-Schrieffer Hamiltonian, and hence corresponds to the 
e — >■ — oo case. The insets magnifies the negative-energy peaks 
of the 7 = 32 fesTx case. The temperature is T = 0.1 Tj^ 

quantify this reduction by studying the appeareance of 
singularities in the resolvents that translate into peaks 
of the spectral function. In order to perform this study, 
we have used the above SCLA applied to the Coqblin- 
SchricfFcr model. 

Let Eq and Ei be the bare energies of the spin zero 
and spin one multiplets. The dressed multiplet energies 
Eq , El are given by solving the equations 

e*^Eo + ^{j:{e*)} 

E*=Ei+di{i:{E*)}, 

where S(cj) is the pseudo-fermion self-energy. Since / is 
equal to the singlet-triplet excitation energy, / = Ei — Eq, 
the threshold shift AI is given by 

AI = (El - E*) - {El -Eo)^^ mE*) - ^(E*)} 

(16) 

Figure ini shows the absolute value of the threshold shift 
as a function of the excitation energy / in units. The 
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FIG. 6: Shift of the excitation energy, |AJ|, versus excitation 
energy 7, crosses. All resuhs are in Tj^ units. F ull line, fit of 
the numerical results to the function |A/| — \/ ax'^ -\-bx -\- c 
where x — I/ln(I/kB T^). Hence, for values of / larger than 
the ones of the present graph, the threshold shift, | AJ| , follows 
the asymptotic behavior ~ I /ln{I /ks 2x). 



function |A/| — \/ax^+bx + c with x — I/ln{I) is an 
excellent fit for a large range of values of /. The fitting 
parameters are a = 0.0029, h = 0.6700 and c = —3.3361. 
Hence, for asymptotically large excitation energies, |A/| 
is of the order of ^/a I /ln{I /ks T^). This result is com- 
pletely equivalent to the asymptotic behavior found for 
the shift of Kondo peaks under a magnetic field B as es- 
timated by Rosch and co-workers^ and by studying the 
shift of spinon density of states under a magnetic field by 
Moore and Wen^. 



IV. DISCUSSION AND CONCLUSIONS 

The above results show that many-body effects in mag- 
netic lETS can be readily studied with usual Kondo- 
physics tools such as the self-consistent approaches 
NCA (for the Anderson model) and SCLA (for the 
Coqblin-Schrieffer one). Comparison with one-electron 
approaches^ is in overall agreement and permits us to 
discern many-body effects that appear as a consequence 
of the impurity excitation's coherence when the inelas- 
tic channels open. Here, we have studied the case of 
a singlet-triplet excitation, when the tunneling electrons 
have enough energy to overcome the singlet-triplet energy 
difference, /. This situation has been achieved exper- 
imentally in carbon nanotubesi^. Related experiments 
are the ones performed in Mn dimers^. However, we ex- 
pect these results to be of relevance for many lETS cases 
when the energy scales correspond to the Kondo ones. 

The one-electron magnetic lETS theory runs parallel 
to the NCA. Namely, a spin excitation takes place via 
a charge transfer process. During this process the mag- 
netic impurity changes its charge state. The decay of the 
charge state back into the adsorbed state can lead to a fi- 
nal state different from the initial one. This is the essence 
of the excitation process. NCA builds on the same idea 



in a self-consistent way such that charge exchange pro- 
cesses are included to all orders, keeping their coherence. 
Hence, Kondo-like features are included. 

The many-body features revealed in this study can 
be summarized by the distorsion of the spectral prop- 
erties of the impurity as compared to the one-electron 
case. Two peaks appear close to the inelastic thresh- 
olds that are due to spin fluctuations when at least two 
spin states become degenerate similarly to the 'usual' 
Kondo case. The inelatic peaks are shifted with respect 
to the energy degeneracy point (the inelastic threshold) . 
This leads to a narrowing of the inelastic gap. This nar- 
rowing increases as / increases following the asymptotic 
law |A/| I/ln{I/kB T^)- This is exactly the behavior 
found for a Kondo peak split by a magnetic field^. 

The resemblance of the present results with those for 
the Kondo effect in the presence of magnetic fields^"— 
is due to the similarities between the two physical pro- 
cesses: in both cases, there is a magnetic excitation, in 
the present case due to the interaction between two local- 
ized spins, and in the magnetic-field case due to Zeeman 
energy splitting, and when the electron energy is large 
enough to open the excited channel, the ground and ex- 
cited states are connected via inelastic spin-fiip electron 
collisions (note that in the present case involving singlet- 
triplet transitions, these are both of spin-fiip and non- 
spin- flip type). 

Finally, the many-body effects described here should 
be observable at large enough Kondo temperatures, T]^. 
This implies that the molecule should be in the Kondo 
regime (|e| < T) but with a sizable T°. (i.e. |e|/r - 0(1) 
). In this situation, peaks at the inelastic thresholds will 
be of many-body nature, leading to a strong renormal- 
ization of the thresholds (see Fig. [S]) which can have im- 
portant consequences in the use of magnetic lETS as a 
spectroscopic tool. The actual observation of many-body 
features may however depend on the measuring proce- 
dure. STM measurements involve several orbitals that 
may or may not be involved in Kondo physics. For real- 
istic systems, the identification of spectral function with 
measured conductance is sometimes not straightforward. 
Although the Kondo peak may prevail in the spectral 
function, the multi-orbital character of the STM conduc- 
tance can lead to channel interference and other effects 
with the consequent appearance of complex Fano pro- 
files^i^. 

Our results show that once Kondo-like features are 
present, they are more robust than usual Kondo peaks. 
Indeed, while Kondo peaks completely disappear at tem- 
peratures a few times the Kondo temperature, inelastic 
many-body features survive in this temperature range, if 
the temperature is smaller than the excitation energy. 
Actually, two energy scales determine the many-body 
properties of lETS: the coherence scale given by ks T^^, 
and the excitation energy. Hence, not surprisingly as 
noticed by Hurley and coworkers^!, some experimental 
lETS show spike-like features at the lETS thresholds, 
where the results of our present work should be taken 
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into consideration. 
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